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Abstract 

We compute next-to-leading order (NLO) perturbative QCD corrections to the 
correlators of interpolating pentaquark currents. We employ modular techniques in 
configuration space which saves us from the onus of having to do loop calculations. 
The modular technique is explained in some detail. We present explicit NLO results 
for several interpolating pentaquark currents that have been written down in the 
literature. Our modular approach is easily adapted to the case of NLO corrections 
to multiquark correlators with an arbitrary number of quarks/antiquarks. 



1 Introduction 

The discovery of exotic quark states and bound states of gluons would be another manifes- 
tation of QCD, allowing for a quantitative check of its features and finding numerical values 
of some important QCD parameters. While glueballs are certainly the most searched for 
states in QCD [1], there is also much interest in exotic states in strong interactions, i.e. 
states built from quarks within QCD which differ from the simplest valence quark content 
of mesons or baryons (see e.g. Refs. [21 [3]). The theoretical investigation of multiquark 
states {Q^Q"^, n + m > 3) and the experimental search for them may provide important 
information on the properties of the interaction of quarks and gluons at large distances. 
Up to recently major efforts have been directed to the study of the dibaryon spectrum 
{n = 6, m = 0), both theoretically and experimentally |4]. This particular six-fold state 
is rather peculiar as it is close to the deuteron, building bridges to applications of QCD 
to medium-energy nuclear physics [5l|6]. In Ref. |2], Jaffe predicted that there might exist 
a stable six-quark S-wave state - a dihyperyon H - which is a singlet with respect to 
both colour and flavour SU{2) (with strangeness —2) with the quantum numbers J^ = 0^ 
and a mass around 2150 MeV. The quantum numbers of the H state are identical to the 
quantum numbers of the (AA) pair of two A(1115) hyperons, and its mass is smaller than 
the sum of the masses of the two A hyperons. The H state is therefore stable with respect 
to strong interactions and can decay only through weak interactions. To the best of our 
knowledge the famous dibaryon state H is the first state to attract attention in the mod- 
ern context of QCD. Thereafter there were efforts to identify some mesons in QCD (scalar 
mesons as a KK molecule) with a four quark state in order to explain their properties and 
in particular their production and decay patterns [7J . In the intervening years the interest 
in exotic states has mainly shifted to tetraquarks and pentaquarks. 

The study of bound states in QCD is a difficult problem. After almost 40 years of 
research it is clear that the most promising approach is very likely given by lattice QCD in 



particular since the computer power and computer algorithms have advanced much since 
the first introduction of lattice QCD in the early seventies of the last century (results are 
given for instance in Ref. [8]). Besides lattice QCD, model dependent approaches have 
been used in Refs. [9l |T0] , for example in the framework of the MIT quark-bag model [11] . 
It is important to test these model predictions solely on the basis of fundamental principles 
of QCD. Such a test can be made by means of the method of QCD sum rules, using either 
the technique of finite-energy sum rules [12] or that of Borel sum rules [13] . 

The operator product expansion and QCD sum rules serve as a solid testing ground 
for many calculations in the theory of hadrons. The method of QCD sum rules is based 
on the fundamental field theoretic principles of QCD, and has proved its effectiveness 
in calculations of the masses of mesons [T3| HH |T8] and baryons [T5t [T6| ITT] . However, 
the reliability of perturbative calculations requires a thorough check in particular in the 
uncharted territory of exotic multiquark states where the collected experimental material 
is rather small. It is therefore worthwhile to compute some examples in order to get a 
feeling for the structure of the perturbative series. Work in this direction is under way. 

Glueballs have been previously analyzed in the context of QCD sum rules. The pertur- 
bative QCD corrections to the sum rules were found to be very large pjj. Exotic mesonic 
states have been analyzed in Refs. [20j. QCD sum rules for ordinary three quark baryon 
states have been widely studied. In particular, the correlators of baryonic currents with fi- 
nite mass heavy quarks have been calculated at next-to- leading order of perturbative QCD, 
allowing for further improvements in the precision of QCD sum rule predictions [21]. It 
is known that NLO perturbative corrections to baryon sum rules are large [18], El]. They 
are expected to be even larger for multiquark states with n > 3 quarks. Different aspects 
of such n > 3 multiquark states in QCD have already been discussed some time ago [22] . 
One feature of n > 3 multiquark states is that they have a large internal weight of colour 
states [3]. 

The immediate purpose of the present investigation is to concentrate on a type of exotic 



multiquark state called pentaquarks - states with baryon quantum numbers that contain 
an additional quark- ant iquark pair. These states have been discovered experimentally 
by different collaborations: LEPS collaboration (Japan) [23], DIANA Collaboration (Rus- 
sia) 121], CLAS Collaboration (USA) [25], and SAPHIR Collaboration (Germany) [26j. The 
results of the present investigation will open the possibility for a high-precision description 
of these experimental data on pentaquarks. The investigation is also important for further 
experimental precision studies on these and related states at DESY (HERMES Collabo- 
ration [27]) and CERN (NA49 Collaboration [28j). There is also a proposal to launch an 
experimental study of pentaquark baryons at meson factories [29j. Experimentally these 
collaborations are using different apparata and techniques but theoretically the observed 
states should be understood within QCD. While the first principle numerical computa- 
tion on the lattice gave rather positive results [30], analytical methods and in particularly 
method of QCD sum rules should definitely be developed for a reliable identification of the 
new states in the hadronic spectrum. 

It appears that the experimental confirmation of these states is problematic at the 
moment as some collaborations have reconsidered their results and conclusions. However, 
there is no doubt that such states are possible within QCD and the theoretical study should 
continue. In case of a definite positive indication from theory the experimental searches 
could certainly proceed in a much more efficient way. 

In particular a dedicated experiment has given a negative result in the direct search of 
the pentaquark state [31]. A review of the present experimental situation can be found in 
Refs. [32] (see also Ref. [33]). Note that more lattice studies have become available [34], 
some with a negative outcome as concerns the existence of pentaquark states [35j. This 
makes the task of the theory even more challenging. Either one has to show that such 
states do not form for some reason, or to suggest a new mass scale of these states and 
to identify the appropiate decay modes for their determination [36]. The first task is 
difficult in as much as one touches on the problem of bound state formation and therefore 



of the (confined) strong coupling. The latter problem ultimately requires the calculation 
of perturbative corrections to the operator product expansion used within QCD sum rules. 

Sum rule calculations of pentaquarks and corresponding critical analysis' have been 
presented in numerous papers [371 EHl EH]. While the accuracy of the QCD sum rule 
method is about ~ 20% at present, the results obtained agree with experimental claims 
and model predictions [lOj. However, within the QCD sum rule method it is not possible 
to predict whether the mass of the lowest pentaquark state lies above or below the Kp 
threshold (i.e., whether it is stable). 

The aim of the present paper is to create a framework for an accurate sum rule analysis 
of the properties of pentaquark states. The study of pentaquark states within the QCD sum 
rule method requires a precise knowledge of the correlators of the pentaquark interpolating 
currents. In this paper we present perturbative next-to-leading order calculations of the 
relevant correlators in QCD. 

2 NLO corrections to the correlation function 

According to the QCD sum rule approach to hadron properties the principal quantity to 
be analyzed is the correlation function 

n(g) = z|d^xe^'^^(0|Tj(a;)j(0)|0), (1) 

where j{x) is a local current operator with the quantum numbers of the hadron state, 
termed the interpolating current of the hadron state. The construction of the conjugate 
operator j{x) depends on whether the hadron is a fermion or a boson. For fermionic states 
such as the ordinary baryon states or the pentaquark states dealt with in this paper one 
has j{x) = j"'^(x)7°. For bosonic states (mesons, tetraquarks, . . . ) the conjugate operator 
is just the adjoint operator, j{x) = j^(x). The result of the sum rule analysis depends 
strongly on the choice of the interpolating current as has been shown already in the case 
of the dibaryon state [22] . 



In the case of a given pentaquark state the pentaquark current j{x) is a local scalar 
current with the quantum numbers of that pentaquark baryon. For instance, take the 
ground state pentaquark state O^. The current is constructed from five quark fields, such 
that its projection onto the real pentaquark baryon state |0'*"(p)) (within the assumption 
that this state exists) is nonzero: 

(0|j(0)|e+(p))=Ae+, p' = ml,. (2) 

Since such a current j{x) is not unique, the question of its optimal choice arises immediately 
(see Appendix B for a discussion of this issue). We recall that the problem of choosing the 
current already arose in the case of baryons fibi [T6] where the currents are constructed 
from three quark fields. When the current is constructed from five quark fields as in our 
case, this problem is much more complicated, since the number of independent currents 
with the given quantum numbers is much larger (see also Ref. tSSj)- The treatment of 
the current j{x) in the most general form, i.e. in the form of a linear combination of all 
the independent local operators with the quantum numbers of 6, is a very cumbersome 
problem. Therefore we confine ourselves to the choice of a few of the simplest currents 
with the required quantum numbers and analyze the dependence of our results on the 
properties of these currents. As in the case of mesons and baryons, we shall construct the 
current j{x) from quark fields without derivatives. 

In accordance with the method of QCD sum rules, we shall calculate the correlation 
function ([T]) by means of Wilson's operator expansion, assuming that the vacuum expecta- 
tion values of the local operators (the so-called condensates) are nonzero. The calculations 
must be performed in the eulidean region — g^ > 1 GeV^. In this region, the effective 
strong interaction constant as is not very small and one has to calculate the coefficient 
functions of the operator expansion at least at NLO in perturbation theory. Therefore, 
the correlator function 11 should be calculated at NLO in a^. The fact that this may be 
necessary for calculations of physical quantities in the framework of the sum rule method is 
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Figure 1: LO contribution (a) and examples for the NLO propagator (b) and dipropagator 
corrections (c) 



confirmed by previous applications of the sum rule method, in particular by the calculation 
of baryon masses. 

The LO calculation falls into the category of the sunset diagrams (cf. Fig. [1^). Sunset 
diagrams are directly calculable in configuration space [IQ]. These types of diagrams also 
appear in the effective low energy gluon correlator for light quarks [41]. The corrections 
are of two types. The propagator-type corrections depicted in Fig. [lb are straightforward 
and can be easily added with no effort at all. The second type of corrections depicted 
in Fig. [It correspond to the irreducible diagrams of the fish type. They are rather well 
known in the massless limit (the more complicated massive case was analyzed in Ref. |42]). 
In order to deal with the diversity of interpolating currents that have been proposed in 
the literature for the pentaquark states, we have developed a modular calculation method 
in configuration space. The modular method reduces the perturbative calculations of the 
present paper to pure algebraic calculations ^3] . 



3 Presentation of the modular method 

As already mentioned, the two required main modules of our method are the propagator 
correction Si{x) and the dipropagator correction 5*2 (x) which read {t = 7^a;^) 

^i(a:)|NLO = 5i(a;)|Lo{l-C^^(/i^x2)'} = ^o(x') {l - C^^ (^^x^)'} , 



(3) 



^2(x)|nlo = SoixyU^t-e^t'^^ifiyyx 



Y^Y 



1 ii\ n u ,2 



e^Yj ^"^^ + 1,^ + 2 ' "^ ^^'^ 



+ (^ + 4) rf ^ ® Ts „Am^. (4) 



where in the Euclidean domain one has 



.2^_ ^r(2-5) 

2vr2-e(a;2^2-e- 



The renormalization scale ^^ is appropriate for calculations in configuration space if one 
wants to avoid the appearance of ln(47r) and 7^; terms. The scale ^x is related to the scale 
Jx of the MS-scheme by 

/i, = /ie^^/2. (6) 

The direct product signs "®" in the dipropagator correction 5*2 (x) serve to distinguish 
between the two fermion lines involved in the gluon exchange. Finally, 

rr = i^'YY^ = liYYY - YYY) (7) 

is the totally antisymmetric product of three gamma matrices. Eqs. ([3]) and (jl]) allow one 
to calculate the corrections to n-quark/antiquark current correlators of any composition 
without having to calculate any integrals. In Ref. |43j we have presented results for a 
model current with five different flavours. In this paper we deal with several interpolating 
currents suggested in the literature including the equal flavour case. Because of flavours 
appearing twice or more times in the interpolating current, the Wick contraction will result 
in a main contribution and different "cross-over" contributions. 
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Before giving our results for the various interpolating currents, we have to deal with 
renormalization. Corresponding to the propagator and dipropagator corrections, the cor- 
relator function is renormalized by the wave function renormalization factor and the vertex 
renormalization factor, respectively. Concerning the vertex renormalization factor one has 
to account for mixing effects. Mixing can occur when gluons are exchanged between quark 
lines in the pentaquark correlation function. Mixing is taken into account through the 
subtraction of corresponding vertex divergences generated by an operator that can admix 
to the initial current. The general formula reads 

{ipi ® ^j)r = i^i ® ^i) - ^{hi' 8) Ijf + -crff ® <,v ) (V^^' ® ^f). (8) 

Here a"'^ = z/2[7",7^] and all numbers are calculated in Feynman gauge. Note that the 
part proportional to a is gauge independent. The renormalization within our modular 
approach follows from the above line of arguments and leads to counter terms which are 
listed in explicit form in the following. 

Once the correlator function is renormalized, we can calculate the spectral density 
corresponding to the correlator. For this purpose, instead of calculating exphcitely via 



niq)=tj d^xe^'^'iOlTjixymiO). (9) 

in momentum space, one can use the formulas given in Appendix A. 

4 Results for pentaquarks of the first kind 

The correlators in this section are 

(0|Tj(a;)j(0)|0) = So{x')\x'fm,{x'). (10) 

We start with different interpolating currents proposed for the lowest pentaquark state 
0+ at 1530 MeV with quantum numbers J^ = 1/2+ and 5 = 1. Ref. [37] gives an 



overview over pentaquarks which are built up by a diquark, a meson and a single quark. 
In the following these currents will be called pentaquark currents of the first kind. The 
interpolating current with isospin J = is given by 
. . 1 



V0[X) = -/l^abc 



u^{x)C-f5db{x) {ue{x)se{x)i-f5dc{x) - {u ^ d)} . (11) 



Due to the two parts of the interpolating current, there are two diagonal and two mixed 
bare contributions, 



Ksi^') - 


I TT Ve 3 / J I TT \e J ) 

-3{i + ^(^^.^)-(-i + f)} = n-,(.^), 


> 


KU^') - 


-«{->^-^^)'G4)}-^{->^-^^K-74)) 


(12) 



The counter terms for the current read 






An- = -18^ --^+^--) = An-. (13) 



The singularities cancel in the renormalized results which reads 
Kr(^') = 180{l + ^(^ + ln(^^x^))}-12{l + ^(| + 31n(^^a:^))} 



n V o 

tt.« /. _. , o ox\ 1 „ r, a, /17 



Mi + ^if^l-'r{^-Hf^y))\ = Kn(-% 



'"noR^^ ^ 



Kr(^') = 18{l + ^(5 + 71n(^^x2))} + 3{n-^(^-ln(^^x^))} = n^ 

(14) 



In order to calculate the spectral density we have to treat the first order correction and 
the counter term separately. The reason is that these two contributions have different 
x^ powers. The general procedure for the calculation of the spectral density is left to 
Appendix A. The result for the spectral density reads 

pis) = ^^^-^ Uo + —(bi + Ci + —Bo + Bo In (^]] ] (15) 

^^ ' 604800(47r)8 \ " tt V 105 V^//J 
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where Aq is the LO contribution, Bq and Bi are the singular and finite parts of the NLO 
result, resp., and Cq ( = — -Bq) and Ci are the singular and finite parts of the counter term. 
Collecting all contributions one obtains 



An 



IT \ e 



vo 



372 + 60— (- + 25 

n \ e 15/ n \ e 



(16) 



The spectral density finally reads 

c5 



p,„(.) = ^^^ (372 + 60^(9 In f^Ul^y 

'^'^"^ ^ 604800 (47r) 8 \ 7i\ \sj 105/ 

31s^ f a. /45 //22\ 7429\1 
50400(47r)8 \ ^ V \31 ^ \7 ) ^ ~6ET J j ' 



(17) 



The perturbative correction can be seen to be rather large, cf. 7429/651 (as/vr). For the 
remaining currents proposed in Ref. [37] one has 



= Tf^'^f"^ 



ri2{x) 



1 

1 
1 

^abc 



^abc 



ul{x)C'j5db{x) {ue{x)se{x)i'y5dc{x) + (u ^ d)} , 
mJ (x)C7''4(x) {7/.75Me(a;)se(x)i754(x) -(mo rf)} 



a6c I Mr(x)C7''M(,(x) -f^-f5de{x)Se{x)i-f5dc{x) + {u ^ d)j 
ul{x)C-f^Ub{x) -f^-f5Ue{x)Se{x)i-f5Uc{x). 



(18) 
(19) 
(20) 
(21) 



Again we only give results for the spectral densities. They read 

o5 r ^, /i /,,2 



Pv[ (s) 
Pr,;(s) 



2100(47r)8 1 TT U Is/ 105/ 



17s^ 



6300(47r)8 



525(47r)8 



, as ^2255 V 
^ vr U08j/' 

TT \ 6 \ s / 360, 



1 + 



«., 



6 I s 



15877' 



525(47r)8 i" ' vr \ Q '" \ s J 10080 
The perturbative corrections can become as large as 2255/408(a;s/7i") 



(22) 
(23) 
(24) 
(25) 
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5 Pentaquarks of the second kind 

Pentaquarks of the second kind consist of two diquarks and one antiquark. Three possible 
choices are given in Ref. [30], and, with some small deviations, in Refs. [HlllS]. They read 

9f{x) = eabc^aefebgh u^{x)Cdf{x) u^ {x)C-f5dh{x) CsJ(x), (26) 

Op^ix) = eahc^aef^bgh ul{x)C'-f^df{x) u^g{x)C'^^''^^dh{x) Cs'^{x), (27) 



et'^ix) 



(^abc^aef^bgh 



u^{x)Cdf{x) u {x)C-f^-f5dh{x) 75Cs^(x). 



^7;Tt 



(2J 



First we consider the case when the Lorentz index /x in the correlator is contracted. One 
then has 






1575(47r)8 
1575(477)8 



[s . 



(29) 
(30) 



The perturbative correction to pg-^ is the largest of all the cases treated in this paper with 
1021/70(0:5/77). When the Lorentz index is not contracted, we obtain an ordinary and a 
cross-over contribution for the correlators 



u^,:{x) = {o\Tel^%x)et'\o)\o), ^ = 2,3. 

which are the same for both currents, namely 

m%"(x) = -38Ax^t{x^g^''-2x^x'')h + — [i2': 

-1 1 



(31) 



'-9,3 



,2 2 



i^- 



+ 256x'rir-{^^lxJ(- + l), 



nf^"(x) = 384x'it(xV"-2x'^x")— f/i^x^) 



1 1 

£ ~ 2 



-192x'YtY^{f^Wy[l + l 



(32) 



(33) 



The total contribution is 

U'^'^six) = -384x2^(xV'-2x''x")|l + — (/i^x^)' 

+Q4x'r:tr^{fiix'y(^-^ + iy 
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8 31 



(34) 



The counter term reads 



An^r(x) = ^- 

vr e 



(y 1 
■' 896x2^6 (a; V - Sx^a;'^) - 64xVi^7'' • (35) 



The evaluation of the absorptive parts of (1M|) and (135|) which are needed to obtain the 
corresponding spectral functions are technically involved and will not be done here. 

6 Conclusion 

The quark condensate (qq) is the only dimensionful quantity number that appears in the 
sum rule analysis since we assume factorization for the vacuum expectation value of the 
six quark operators |l6]). Inclusion of terms ~ rus should not substantially change the 
quantitative results as the mass of the strange quark is small [47j. It therefore follows 
from dimensional arguments that the mass of the pentaquark baryon is mg ~ i\{QQ)\Y 
as long as power corrections determine the mass. Consequently, me should not change 
by more than ~ 10% if (qq) varies by 30%. Such a variation is quite possible because of 
the uncertainties in the light quark masses determined from the numerical value of the 
light quark condensate as calculated from the PCAC relation for the pion. The analogous 
expression for the strange quark condensate obtains some corrections due to the s-quark 
mass which are well under control |48]. Nevertheless, this still leaves the uncertainty 
whether the pentaquark state is above or below the threshold. 

We recall in this respect that the high accuracy of the MIT quark-bag model permitted 
Jaffe to conclude that the dibaryon state H lies below the AA threshold and is therefore 
stable. The same conclusion was drawn from a model calculation based on chiral solitons 
in Ref. [10]. However, for a model independent approach, the relatively low accuracy of 
the method in the determination of the mass (~ 15%) does not make it possible to draw 
any conclusion about whether the mass of the exotic baryon lies below or above the KN 
threshold. 



13 



In this paper we have calculated NLO perturbative corrections to the correlator of var- 
ious pentaquark currents. We have shown that such a calculation can be done by purely 
algebraic means for any given interpolating current using the modular methods developed 
by us in detail. As it turns out, the NLO corrections to the correlators are large. As the 
coupling constant is large at the relevant energy scale |19], the large perturbative a^ correc- 
tions will heavily change the relative weight of the perturbative and the nonperturbative 
condensate terms. Such large NLO corrections to the current correlators could possibly 
introduce an element of uncertainty and unreliability in the QCD sum rule analysis of 
pentaquark states. 
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A Explicit derivation of the spectral density 

The momentum space representation of the correlator function can be obtained from the 
configuration space representation by using the integration formula 

U{p) = 27r^+i r (^] Jx{px)U{x)x^^+^dx (Al) 

where X = 1 — e and J\{z) is the Bessel function of the first kind. In the case that the 
correlation function Il{x) = (x^)^" is a simple power in x^, the integral can be obtained 
analytically. The result is 



n.(/) = n'^' ( ^ V ^^^';/, ^^ (A2) 



4 i Via) 
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The corresponding spectral density is given by the discontinuity of the correlation function 
where the discontinuity of the correlation function (in the Euclidean domain!) lies along 
the negative real axis. One obtains 

In order to calculate the spectral density, we have to use the scalar correlation function. The 
vector-type correlation function of the pentaquarks (as well as those of all states composed 
of fermions) can be obtained from the derivative of this scalar correlation function F(x^). 
One has 



dx"^ 



d,F{x') = 2x,^^^ = xj{x'). (A4) 



Given the function /(x^), the scalar correlation function is obtained by integrating over 
x^/2. In case of pentaquarks, we have 

fix') = {So{x')y {x^f [a + ^{^^lxyB] (A5) 

where A contains the LO contribution and the counter term and where B contains the 
NLO contribution. Recalling the x^-dependence of 5*0(0;^) in Eq. ([S]), one obtains 

F{x^) = \{^^i^) j{{x'r''A + ^[^^t)\x-r-'B]dx- = 

1 /-r(2-£:)\ (, o,.,_7 , a.,/ n\e, 9^«=._Q„5£: - 7 



[{x'r-'A+^{^j?:)\x-r-'B-^ 



2{he-l)\ 27r2-£ / ^ -^ ^^^ 65 - 7 



2\\^ 



iSoix')) 



2(5£ - 7) 

For the corresponding spectral density one has 
,, TT^-- f -n2-e) Y 



{A + ^{,yyB[l^^)}. (A6) 



(s/Ar-^^A , as ifiins/Af-'^B f^^e 



r(7 - 5£)r(6 - 4£) TT r(7-6£)r(6-5£) V 7 
-7r2-^(s/4)5-4- /-r(2-£)\^ 



2r(8 - be)V{Q - As) \ 27r2-£ ) 
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X 



A + ^l^] B (l + i) (l + (*.(7) + ii6) + 27e) e 



2r(8 - 5£)r(6 - Ae) \ 2n^-^ J { n \ s J V 7/ V 15 

(A7) 

where we have made use of the expansion r(a — e) = T{a) (1 — eipla) + 0(6^)) and where 
we have incorporated the scale change /i^. = e'^^/2/2. Here ■?/'(a) = r'(a)/r(a) is the 
polygamma function. We then use 

a,, /Co \ n -Bo 



A = Ao + — [ — + Ci), B = ^ + Bi (A8) 

where Aq is the LO contribution, Bq and i?i are the singular resp. finite NLO contribution, 

and Co and Ci are the singular resp. finite contribution of the counter term (Co = — -Bo)- 

One finally obtains 

^7r^-^(g/4)5-4^ / -r(2 -£)^' 
2r(8 - 5£)r(6 - Ae) I ~2^ 



Pf{s) = .T./0 .^^T^^a — 1 o o_e >< 



^Ao + ^(B, + C^ + ^Bo + Boln('^]]]. (A9) 



604800(47r)8 [ vr V 105 \s 

Because the singularities cancel one can set £ = in the last step. 

B On the choice for the interpolating current 

Formally, the QCD sum rule method dictates a priori the only condition for the choice of 
the current: it must possess the required quantum numbers. However, a posteriori (after 
the fit) positive results can be obtained only for "physical" currents. In particular, the 
two requirements formulated below are usually necessary. This may be an indication that 
there exists a criterion which makes it possible to select the optimal ("physical") current 
from the set of currents with the quantum numbers of the given channel. This argument 
was also given in Ref . [5] , in which the deuteron mass was calculated by the same method. 
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The choice of the interpolating current is crucial in this respect and has to be considered 
very carefully. The "physical" current G satisfies the following two requirements. First, 
there exists a nonzero nonrelativistic limit for it (i.e. if the quark field iIj{x) is represented 
in the standard manner in terms of the small and large components, the term containing 
only the large component will be nonzero). We note that in Refs. [151 EO] it was already 
pointed out that the existence of a nonrelativistic limit is a desirable property for the 
construction of currents when employing the QCD sum rule method. The demand for the 
existence of such a limit is quite natural as the results should be reproducible in some 
effective potential model of constituent quarks. Second, the colour (and flavour) structure 
is important. This will be explained in some detail in the following. 

According to Ref. [30], "the qqqqq state can be decomposed into a pair of color singlet 
states as qqq and qq [■■■]■ For instance, one can start a study with a simple minded local 
operator for the G'^(1540), which is constructed from the product of a neutron operator and 
a K~^ operator such as 6 = eabddaClbUbjdc^Se'j^Ue). The two-point correlation function 
composed of this operator, in general, couples not only to the 6 state (single hadron) but 
also to the two-hadron states such as an interacting KN system. Even worse, when the 
mass of the qqqqq state is higher than the threshold of the hadronic two-body system, 
the two-point function should be dominated by the two-hadron states. Thus, a specific 
operator with as little overlap with the hadronic two-body states as possible is desired in 
order to identify the signal of the pentaquark states [. . . ]." And following Ref. [37], "isospin 
and color structure guarantee that these currents will never couple to a K~^n molecule or 
any other K~^n intermediate state [...]." This is the reason to use a non-trivial colour 
structure in the previous paper. 

We briefly comment also on a second choice of currents with non-trivial flavour structure 
which was not considered in this paper. If we represent the current B in the form of a 
product "singlet®singlet" with respect to colour, i.e., in the form B = \E'(a:)r$(x), where 
each of the operators \1', $ is a colour singlet and F is a string of Dirac 7 matrices, we can 
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choose 

where the colour-singlet operator \E' is a flavour octet 



Q = ^it{xMl{x) (Bl) 



Here a is a spinor index. The flavour octet ( lB2p has the quantum numbers of the baryon 
octet and has e.g. been used in Refs. [151 US] to calculate the properties of light baryons 
within the QCD sum rule method. Thus, if the current B is represented in the form 
"singlet®singlet" with respect to colour, it has to have the structure "octet®octet" with 
respect to flavour (each colour singlet is a flavour octet). Physically, the following picture 
emerges. The colourless state 9 splits into two colourless clusters which separate at large 
distances to become a meson and a baryon. As a result, we conclude that the current G 
is the most "physical" current in the sense that it has a nonrelativistic limit and that it is 
constructed as a "octet®octet" state with respect to flavour. 
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